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Abstract 

The DFT representation is applied, to a ordered-chaotic finite-duration se- 
quences set generated by the Henon mapping, for describing the order-chaos 
transition. This representation has the advantage that it may be applied to 
a relatively shorter discret-time series. The bifurcation diagram is produced 
and the largest Lyapounov exponent is calculated for each time series of the set, 
showing a good agreement with the results obtained by the DFT representation. 
The threshold value of the control parameter for the order-chaos transition was 
determined to be A=1.0532. 



1 Introduction 

The Discrete Fourier Transform (DFT) representation method has been yet ap- 
plied for describing the order-chaos transition in a dynamical model represented 
by the Poincare mapping ||. This representation has the advantage that can 
operate over a relatively shorter discrete-time series, and it was demonstrated 
that the DFT representation of a chaotic discrete-time series is possible, never- 
theless this discrete-time series type is to be a nonlinear one. Literature reports 
that transition from a regular to a chaotic bahavior has been studied in many 
dynamical system, including bilogical systems [^[ || . 

In this work, the Henon mapping dynamics is considered, and the DFT 
representation is applied to describing the order-chaos transition. The Henon 
mapping is a dynamical model widely used that can produce a chaotic discrete- 
time series, and therefore it is important to evaluate the largest Lyapounov 
exponent for each of the discrete-time series produced. 

In 1968, Michel Henon || proposed a quadratic mapping of the plane as 
a way for stuying a dynamical system such as the motion of a asteroid or a 
charged particle in an accelerator. Henon mapping is based on the idea that 
one can reduce the study of a conservative system with two degree of freedom to 
the study of area-preserving mappings of the plane. Henon mappings simulate 
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the behavior of a physical system and show that many such systems are more 
complex than previously imaged. The nature of this complexity just begins to 
be understood now a day by researchers. 

Wolf et al |Hj| included the Henon mapping Qin their study for determining 
Lyapounov exponents from a discrete-time series. The Lyapounov exponents 
spectrum was calculated to be Ai = 0.603 bit/iter, and A 2 = —2.34 bit/iter. 
The Lyapounov dimension was found to be Dj = 1.26. 

A control technique has been developed by Ott et al ||. It is based on 
the fact that there are an infinite number of unstable periodic orbits embedded 
within a chaotic attractor. This method is a powerful experimental tool that 
models the dynamics of a system directly from its discrete-time series. The Ott- 
Grebogi-Yorke technique has been used to control many different chaotic system. 
The principles of this technique have been used by an another technique called 
tracking to stabilize underlying unstable periodic orbits in nonchaotic systems. 
To investigate the feasibility of using this method, Christini and Collins j| 
considered the Henon mapping with additive noise, as given by the expression 
x[n] = 1.0 — A * x[n — l] 2 + B * x[n — 2] + £[n], where A is an adjustable 
parameter, B=0.3, and £[n], is a Gaussian white noise with zero mean and 
standard deviation known. They used A as a control parameter, varying into 
the interval [0.2,1.2], to produce a bifurcation diagram of the noise-free model. 
It was their interest to kepp the Henon mapping out of the chaotic region by 
setting A to 1.00. At this condition, the Hcnoin mapping exhibits a stable 
period-4 cycle. 

Khadra et al |^| applied a statistical approach to chaos identification in 
discrete-time series. The method was applied to numerical data generated by 
chaotic systems and to heart rate variability signals of normal subjets and heart 
transplant recipients. To demonstrate the power of this technique they applied 
it to random and chaotic data. The random data were generated from an AR 
random process defined by x[n + 1] = 0.9 * x[n] + y[n], where y[n] are white 
Gaussian random numbers with zero mean and unit variance. The chaotic data 
were generated from the Henon mapping defined as above, using A=1.4, and 
B=0.3. 

Christini et al H developed a real-time, model-independent control tech- 
nique for chaotic and nonchaotic systems. For testing this technique they 
considered the Henon mapping with an adjustable parameter A n , as given by 
x[n] = 1.0 -A n * x[n - l] 2 + B * x[n - 2~l where B=0.3. 

In a recent work of Iasemidis et al |7| about epilepsy, the voltages of the 
recorded signal x(t) at several points in time [t, t + r, • • • , t + (D e — 1) * r], sep- 
arated by a given time delay r were converted to a D e — dimensional vector 
X(t). The process was repeated by increasing t to obtain successive points in the 
phase space, creating a D e — dimensional structure. The phase space portrait 
can be analyzed for determining the presence of an attractor and its dimension. 
The presence of an attractor of low dimension means that of the many modes of 
the system only a portion of them are active during that period of observation. 
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Thus, the effective phase space portrait dimension is much smaller than the full 
dimension of the equations governing the dynamics of the system. They found 
that the epileptogenic focus generates signals characteristic of a nonlinear dy- 
namical system, and that analysis of signals from the centre of the epileptogenic 
focus demonstrated the presence of strange attractor of low dimensionality and 
with positive value for largest Lyapounov exponent (LLE) during the inmediate 
preictal period, during the seizure, and following the postictal period. They 
used the phase-space portrait method for each stage of the seizure process, and 
calculated for each stage the LLE of the signal. 

The electrocorticogram (ECoG) is a highly nonstationary signal. This con- 
dition must be taken into account by the authors in order to modify the Wolf 
algorithm of the LLE estimating. In this frame, The Henon mapping was used 
for testing the LLE estimating algorithm. They considered the Henon mapping 
as given by 



with A=1.4, B=0.3, and the initial condition to be (x[l],j/[l]) = (0.631,0.189). 
They chosen N=2000 data points, with D e —2, and r=l point. The LLE esti- 
mated value was 0.577 bit/iter, that is close to that obtained by Wolf et al, that 
is, 0.603 bit/iter. 

Dynamical control of excitable biological systems is often complicated by 
the difficult and unreliable task of precontrol identification of unstable peri- 
odic orbits. However, chaos control techniques have been applied to a number 
of excitable biological systems. One particular control technique, proportional 
perturbation feedback control, uses isolated electrical stimuli to cause the cells 
to fire at a specified time, thus directly altering the variable of interest, in 
interexcitation interval. Christini et al M illustrated the stable manifold place- 
ment characterization and tracking technique using the chaotic Henon mapping 
defined as given by x[n + 1] = A + B*x[n — 1] — x[?i] 2 , where A—1A and B=0.3, 
and x[n] representing the nth interexcitation interval. Here, the chaotic Henon 
mapping was chosen because its uncontrolled dynamics have no stable periodic 
orbit and because its unstable flip saddle unstable periodic orbit is of the same 
type as that reported experimentally in studies of excitable media spike trains. 

Literature reports many other studies that include the Henon mapping model 
to produce a chaotic time series. Generally, the control parameter value chosen 
for producing a chaotic discrete-time series is A=1.4, but another value that 
satisfies the condition 1.0532< A<1.4 may be used, for which all them the 
corresponding LLE is positive. 

The aim of this research is to evaluate the ability of DFT representation 
of a discrete-time series for describing the order-chaos transition in the Henon 
mapping dynamics, where the control parameter A can be adjusted in order to 
obtain a ordered or a chaotic regime as one decides, and comparing this method 



x[n + 1] = 1.0 - A * x[n] 2 + B * y[n], 
y[n + 1] = B * x[n] 



(1) 
(2) 
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with the LLE method. 

In this work, a 120 discrete-time series set, with 1024 data points in each 
discrete-time series, was produced. For each discrete-time series, its 1024 DFT 
representation was obtained. By observating either the amplitude spectrum and 
the phase spectrum behaviors with the control parameter A of the system, a 
well defined value of this control parameter may be determined, beyond which 
the behaviour either the amplitude spectrum and the phase spectrum becomes 
quite irregularly, in contrast with the regular behavior of these magnitudes in the 
region before the indicated value above. There exists a complete correspondence 
when this result is compared with that obtained from the bifurcation diagram 
and from the LLE plot, both as a function of the control parameter A. The 
threshold value for order-chaos transition for the control parameter A is 1.0532, 
which corresponding to a LLE value close to zero. 

2 Methods and Materials 

2.1 Obtention of the discrete-time series set 

A Henon mapping [[| is an area-preserving map of the plane given by 

x[n + 1] = x[n] * Coscf) — (y[n] — x[n} 2 ) * Sirup, (3) 
x[n + 1] = x[n] * Sirup + (y[n] — x[n} 2 ) * Cooscf) (4) 

Here 4> is the phase angle, chosen as a fixed constant on the interval [0,7r]. An 
initial point (x[l],y[l]) is chosen each times to generate one orbit of the mapping 
by applying a iterative process, and modifying the initial point until a typical 
mapping containing 15 or 20 orbits is produced. It is important to choose an 
adecuate starting point for generating the orbit. 

Using the equations above, a typical-38 orbit Henon mapping was obtained. 
For generating orbits, two linear space for the starting point was created, that 
is, 

x[l] = linspace(0.098, 0.718, 38) (5) 
y[l] = linspace{0mi, 0.581, 38) (6) 

The phase angle </> was chosen as 0=1.111 rad. A starting point is chosen 
each times to generate a 1024 point-orbit of the mapping. 
The Henon mapping can be too formulated as given by 

x[n + 1} = 1.0 - A* x[n} 2 + B *y[n], (7) 
y[n + 1] = B *x[n] (8) 

where B=0.3, and A is a adjustable control parameter. With the initial 
condition (x[l],y[l]) specified, a discrete-time series with a given length by a 
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recursive action is obtained, each times. The nature of the behavior of a data 
sequence obtained can be modified by chosen conveniently a value of the control 
parameter A. The interval [0.2,1.4] was chosen, and a set of evenly spaced values 
of A with a step (5=0.01 was produced to obtain a family of 120 discrete-time 
series, each of them containing N = 2 10 samples. The building up of the 
histogram for each discrete-time series allow us to have inmediatly a simple 
statistical characterization of the discrete-time series. 

In order to make much easier further computations, data points were orga- 
nized into a rectangular matrix containing 120 columns, being each of them a 
time series with the length A = 1024 data points. 

2.2 Determining of the largest Lyapounov exponent Ai 

A system containing one or more positive exponents is to be defined as chaotic. 
For calculating the Lyapounov exponents stectrum some algorithms have been 
developed. One of the most used algorithm that reported by Wolf et al |p"3f . 
From the exponents spectrum, the largest exponent, Ai, decides the behaviour 
of the dynamical system. For calculating the largest exponent, one can too refer 
the algorithm proposed by Rosenstein et al Jll[] . Details of these algorithm can 
be found in refered papers. 

The largest Lyapounov exponent, Ai, for each time series in the set was de- 
termining using a professional software JlJ] . The selection of the time delay, r, 
and the embedding dimension, Z) e , required for reconstructing the phase space 
is part of the problem. The value r=l was chosen considering the optimal filling 
of the phase space method Q . Figure 4 depicts the effect of optimal filling of the 
phase space of the attractor for the value of r indicated above. For embedding 
dimension the value D e —i was chosen. 

2.3 Application of DFT representation: Amplitude spec- 
trum and phase spectrum 

The DFT representation of a finite- length discrete-time series x[n}; n — 0, 1, 2, . . . , A— 
1, is given by 

N-l 

X W = n X[k]e j(2n/N)kn ; n = 0, 1, 2, . . . , N - 1 (9) 

N-l 

X[k] = x[n]e-^ /N)kn ; k = 0, 1, 2, . . . , N - 1 (10) 

n=-0 

The sequence X[k]; k = 0,1, 2,..., A — 1 can be efficiently computed by 
a digital algorithm known as Fast Fourier Transform (FFT), from which the 
amplitude spectrum, |A[fc]| or abs(X[k]), and the phase spectrum, arg(X[k\) or 
angle(X[k]), are both obtained. 
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Next we obtain the DFT representation for each time series in the set by 
the 1024 points-FFT of the set matrix. Since each time series produced by the 
Henon mapping (with A £ [0.2,1.4]) is a finite-duration bounded sequence, it 
satisfies the existence requirement for the DFT representation, that is, 



We considered that is convenient, for some proposes, eliminating the dc 
component of signals produced by the Henon mapping, before the amplitude 
spectrum mesh plot was made. 

2.4 Results 

Figures, from 1 to 3, depict a sample of three representative time series and their 
corresponding histograms. In order to give a better view of each time series it 
was only considered 128 data points of each discrete-time series for plotting but 
it was completely considered when its corresponding histogram was made. In 
either case, the particular value of A parameter was specified. 

The two-dimensional plot of the rectangular matrix produces a figure known 
as Feigenbaum or bifurcation diagram. It reflects the transition of a dynamical 
system towards chaos by a period-doubling process. This can be observed by 
running the program at appendix. Figure for bifurcation diagram is obtaining 
by plotting A parameter values on the horizontal axis and sequence data values 
in the vertical axis. A threshold value A = 1.0532 is defined, beyong which 
the transition of the dynamical system to the chaotic behavior can be observed. 
The bifurcation diagram serves as a reference to control results, and it is usually 
applied in this research frame fijl [ll| . 

Figure 5 depicts results for the largest Lyapounov exponent estimate Ai, 
calculated for each of the time series in the set, as a function of the control 
parameter A. This reference is taken as a control for the discussion of the 
results obtained by applying the DFT representation, and it is too a particular 
aim of this work. 

Amplitude spectrum |A[fc] or abs(X[k]) and phase spectrum arg(X[k]) or 
angle(X[k]), both for the discrete-time series matrix, may be plotted (see ap- 
pendix). In either case, the x axis corresponds to the control parameter A, and 
the y axis corresponds to the frequency, from zero to the Nyquist frequency 
value. In order to obtain the best view of the DFT representation amplitude 
spectrum mesh plot, the dc component on the signal was eliminated. Figure 6 
shows the phase spectrum mesh plot of the DFT representation of the discrete- 
time series set. 
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3 Discussion 



The application of the DFT representation of a chaotic-type signal gives satis- 
factory results for determining the order-chaos transition in the Hcnon mapping 
dynamics. This conclusion is confirmed by comparing that phase spectrum mesh 
plot depicted in Figure 6 with that shown in Figure 5, where the behavior of 
the largest Lyapounov exponent Ai with the control parameter A is shown. A 
critical value for control parameter A can be defined, beyong which sequences 
produced by the Henon mapping exhibit a chaotic behavior, in general. This 
threshold value is close to 1.0532, as is obtained from Figure 5. From the Figure 
5, it can deduce the value 0.618 bit/iter for LLE, which correponding to a value 
1.4 for control parameter A. It is in agree with those values obtained by Wolf 
et al 0.603 bit/iter, and by Iasemidis et al 0, 0.577 bit/iter. 

Some intermediate ordered region in the chaotic behavioral region underly, 
but some authors Jl3| refers them as a consequence of the computational system 
limitations rather than a real behavior of a dynamical system. However, this 
relationship between order and chaos can be observed from the dependence of 
the largest Lyapounov exponent Ai with control parameter, A, whenever the 
largest Lyapounov exponent value tends to zero at bifurcation points embedded 
in chaotic region. The behavior of a dynamical system, for control parameter 
values correponding to largest Lyapounov exponent close to zero, is not chaotic. 
The bifurcation diagram was builded up several times, changing the initial con- 
dition, and changing the step value for the control parameter, and the feature 
of the bifurcation diagram remains. It can be remarked that the order-chaos 
transition for the value A = 1.0532 can be observed either in DFT representa- 
tion amplitude spectrum mesh plot and in DFT representation phase spectrum 
mesh plot. 

On the other hand, Figures from 1 to 3 show the time domain representation 
of some of time series and their corresponding histogram. From Figure 1(A = 
0.99), it can be deduced the bifurcation effect from their histogram. Note that 
for Figures 2 (A = 1.17) and 3(A = 1.24) the control parameter satisfies the 
condition A > 1.0531 in either case, and that Figure 3 corresponds to a value 
for the largest Lyapounov exponent close to zero, for which the discrete-time 
series is not chaotic. It can be remarked that the bifurcation diagram analysis 
is used in studies of the behavior of a dynamical system that tends to chaos 
through period-doubling procedure . 

Applying DFT representation to a discrete-time series in a set is recom- 
mended for qualitative detecting the nature (ordered or chaotic) of the behavior 
of the system, and this can be added to the metric tools of the nonlinear dynam- 
ical analysis methods. Any method that allow us to evaluate the behavior of a 
dynamical system, and that can be applied to a relatively short discrete-time 
series, must be taken into account as a practical method. 
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4 Conclusions 



The application of the discrete-time series DFT representation can be made 
over a chaotic-type discrete-time series, nevertheless this time series is to be 
nonlinear one. This method have neither no hard requirement about the length 
of the experimental data sequence. If one deals with a nonlinear dynamical 
system that can modify its behaviour between order and chaos, a discrete-time 
series produced by the system reflects this change, and the DFT representation 
exhibits in its amplitude spectrum and in its phase spectrum the change in the 
system too. Any method, that allow us to give some criteria about the behaviour 
of a dynamical system and that can operate over a relatively short discrete-time 
series, acquires practical importance. The DFT representation can be used for 
complementing the results obtained by applying the most powerfull methods of 
the nonlinear dynamical analysis. 

5 Appendix 

In order to see the bifurcation diagram, a finest version of all figures including 
the DFT representation phase spectrum mesh plot and the DFT representation 
amplitude mesh plot, the following program may be runned on MatLab. 
(For running on MatLab, write a percent symbol at the begining of each com- 
ment line) 

Comment: To producing a family consisting of 120 time series, each of one 
having 1024 data points. 
Comment: Figures for DFTHM paper 

Comment: To producing a family consisting of 120 time series, each of one hav- 
ing 1024 data points, 
clear 
close all 

Comment: Data 

alpha = 0.2; del = 0.01; nalpha = 121; x = 0.098; y = 0.068; N = 1024; 

beta = 0.3; 

for j = 1 : nalpha; 

Alpha(j) = alpha; 

alpha = alpha + del; 

end; 

for j = 1 : nalpha; 
for i = 1 : TV; 
A(i,j) = x; z = x; 

B{i,j) = y; x = 1 - Alpha(j) * x 2 + y; 

y = beta * z; 

end; 

end; 
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Comment: Figures selected 

j = [80 98 105]; 

for k = 1 : 3; 

alfa(k) = Alpha(j(k)); 

end; 

dispfalpha values are:'); alfa, 

dispfvalues of j equals to 80, 98, and 105 were selected'); 
Comment: Alpha = [0.99 1.17 1.24]; 
pause(l); figure(l); subplot(211); plot(A(l:128,80)); 
title('Figure 1: Time series and its histogram, for A = 0.99'); 
xlabcl('Discrete-time index'); ylabel('Scqucncc values'); 
subplot(212); hist(A(l:128,80),40); xlabel('Sequence values'); 
ylabel('Relative frequency'); 

pause(l); figure(2); subplot(211); plot(A(l:128,98)); 
title('Figure 2: Time series and its histogram, for A = 1.17'); 
xlabel('Discrete-time index'); ylabel('Sequencc values'); 
subplot(212); hist(A(l:128,98),40); xlabel('Sequence values'); 
ylabel('Relative frequency'); 

pause(l); figure(3); subplot(211); plot(A(l:128,105)); 
title('Figurc 3: Time series and its histogram, for A = 1.24'); 
xlabel('Discrete-time index'); ylabel('Scquence values'); 
subplot(212); hist(A(l:128,105),40); xlabcl('Sequence values'); 
ylabel('Relative frequency'); 
Comment:For optimal filling criterium 
pause(l); figure(4); 

title('Figure 4: Optimal filling of the phase space'); 

subplot(221); plot(A(l:N,114),B(l:N,114)); 

titlc('Figure 4a : A = 1.33'); 

ylabel('First difference'); 

subplot(222); plot(A(l:N,116),B(l:N,116)); 

title('Figurc 4b : A = 1.35'); 

ylabel('First difference'); 

subplot(223); plot(A(l:N,118),B(l:N,118)); 

title('Figure 4c : A = 1.37'); 

xlabel('Coordinate'); 

ylabel('First difference'); 

subplot(224); plot(A(l:N,120),B(l:N,120)); 

title('Figure 4d : A = 1.39'); 

xlabel('Coordinate'); 

ylabel('First difference'); 

Comment:Part of the program to obtain the bifurcation diagram. 
pause(l); figure(8); 

for i = 1 : N; hold on; plot(Alpha,A(i,:)); end; grid on; 
titlc('Figure 8: Bifurcation diagram'); xlabel('Control parameter'); 
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ylabel('Sequence values'); 

Comment:Part of the program for largest Lyapounov exponent ploting. 
xr = 0.2 : 0.01 : 0.2 + 120* .01; 



yexl = [ 



-0.309 - 


0.286 - 


0.201 - 


0.219 




0.217 - 0.199 - 


0.180 


- 0.162 - 0.152 - 


0.135. 


-0.111 




- 0.098 




-0.077 


- 0.066 


- 0.044 - 0.032 


-0.013 0.001 -0.034 - 


0.075. 


-0.113 - 





153 - 





196 - 


0.243 




0.294 -0.352 - 


0.419 


- 0.499 - 0.591 - 


0.789. 


-1.002 - 





496 - 





769 - 


1.021 




0.969 - 0.889 - 


0.887 


- 0.827 - 0.897 - 


0.937. 


-0.903 - 





854 - 





838 - 


0.842 




0.865 - 0.862 - 


0.824 


- 0.995 - 0.924 - 


1.013. 


-0.867 - 





634 - 





874 - 


0.917 




0.904 - 1.187 - 


1.233 


- 1.183 - 0.957 - 


0.727. 


-0.602 - 





499 - 





419 - 


0.353 




0.295 -0.243 - 


0.196 


-0.153 -0.113 - 


0.076. 


-0.041 - 





007 - 





043 - 


0.132 




0.241 -0.402 - 


0.642 


- 1.661 - 0.711 - 


0.395. 







232 - 





124 - 


0.041 




0.062 -0.407 - 


0.019 


0.041 0.141 0.188 


0.259. 










0.272 


0.188 





340 0.336 0.363 


0.383 


0.414 0.420 0.455 


0.414. 






0.428 





475 0.447 - 





415 - 0.311 - 0.083 - 


0.006 0.113 0.358 


0.444. 










-0.248 


0.294 





453 0.513 0.499 


0.526 


0.501 0.502 0.549 


0.570. 



0.618] 



pause(l); figure(5); 

plot(xr,yexl,'kx'); hold on; plot(xr,yexl,'k:'); grid on; 
xlabel('Control parameter'); ylabel('Largest Lyapounov exponent'); 
title('Figure 5: Largest Lyapounov exponent behaviour with A'); 
legend('Experimental value','Behaviour',0); zoom on; 

Comment:Part of the program for calculating and ploting the amplitude spec- 
trum mesh and the phase spectrum mesh of the DFT representation. 
A = A — ones(N, 1) * mean(A); k = 9; / = linspace(0,pi, iV/8); 
for k = : N/8 -l;m = 2*fc + l; Ai(k + 1, :) = A(m, :); end; 
Ai = Ai — ones(N/8, l)*mean(Ai); Alpha — linspace(0.2, 1.40, (nalpha— 1)/4); 
for j = : (nalpha — l)/4; p = 4* j + 1; Aii(:,j + 1) = Ai(:,p); end; 
Aii = Aii — ones(N/8, 1) * mean(Aii); ftAii = fft(Aii); 
amftAii = abs(ftAii(:, :)); anftAii = angle(ftAii(:, :)); 
pause(l); figure(7); mesh(Alpha,f,amftAii(l:128,l:30)); 
title( 'Figure 7: Amplitude spectrum mesh'); 
xlabel('Control parameter'); ylabel('Frequency'); 
pause(l); figure(6); mesh(Alpha,f,anftAii(l:128,l:30)); 
title('Figure 6: Phase spectrum mesh'); 
xlabel('Control parameter'); ylabel('Frequency'); 
view(-33,50); 
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Figure 1: Time series and its histogram, for A=0.99 



Figure 2: Time series and its histogram, for A=1.17 



Figure 3: Time series and its histogram, for A=1.24 



Figure 4: Optimal filling of the phase space criterium 



Figure 5: Behavior of the largest Lyapounov exponent 



Figure 6: Phase spectrum mesh 
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Figure 2: Time series and its histogram, for A = 1.17 
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Figure 5: Largest Lyapounov exponent behaviour with A 
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